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ABSTRACT 

We have carried out a systematic study of the properties of the kilo-Hertz quasi- 
periodic oscillations (QPO) observed in the X-ray emission of the neutron star low- 
^ I mass X-ray binary 4U1608-52, using archival data obtained with the Rossi X-ray 

Timing Explorer. We have investigated the quality factor, Q, of the oscillations (de- 
fined as the ratio, ly/ Ai', of the frequency ly of the QPO peak to its full width at half 
maximum Aiy). In order to minimise the effect of long-term frequency drifts, power 
spectra were computed over the shortest times permitted by the data statistics. We 
[ show that the high Q of ~ 200 reported by Berger et al (1996) for the lower frequency 

. kilo-Hz QPO in one of their observations is by no means exceptional, as we observe 

a mean Q value in excess of 150 in 14 out of the 21 observations analysed and Q can 
remain above 200 for thousands of seconds. The frequency of the QPO varies over the 
wide range 560-890 Hz and we find a systematic trend for the coherence time of the 
O ' QPO, estimated as t = Q/[t:u) = 1/[t:Av), to increase with t^, up to a maximum 

|H , level at ~ 800 Hz, beyond which it appears to decrease, at frequencies where the QPO 

^2 ' weakens. There is a more complex relationship between r and the QPO root mean 

squared amplitude (RMS), in which positive and negative correlations can be found. 
A higher-frequency QPO, revealed by correcting for the frequency drift of the 560- 
890 Hz one, has a much lower Q (~ 10) which does not follow the same pattern. We 
. ^ discuss these results in the framework of competing QPO models and show that those 

^ . involving clumps orbiting within or above the accretion disk are ruled out. 

1 INTRODUCTION the QPO discovery paper of iBerger et~al] (|l99(jl . a quality 

factor Q — iz/Au of up to ~ 200 was reported for the QPO, 

Fourier analysis reveals the X-ray emission of low-mass X- together with apparent correlations between its frequency 

ray binaries (LMXBs) to be variable on short timescales, and both its amplitude and width, 
often with several characteristic quasi-periodic oscillation 

1 QPO) peaks ident i fied b e tween ~ 1 Hz and ~ 1 kHz 
van der Kli j Il989l l2000l. |20o3). The highest-frequency 
QPOs, reaching up to ~ 1300 Hz and usually occurring 
in pairs, are of particular interest, as their periods corre- 
spond to the dynamical time-scale in the inner accretion We have retrieved the 1996 and 1998 data of 4U1608-52 from 
disk, where strong-field effects of gravity are crucially im- th e RXTE arch i ve. Th ese d ata have already been di scusse d 
portant. Ther e is no agreement as to the physical origin in iBerger et alJ J1996D and iMendez et all Jl998allbl ll99E 



of the QPOs l|van der Klisll2004h . In the wealth of obser- The 1996 data were recorded during the decline of an X-ray 

vational papers published on QPOs, many have studied the outburst, whereas the 1998 data sampled a whole outburst 

QPO properties as a function of count rate, luminosity, spec- phase from the source. Here, we have restricted the analysis 

tra, lower-frequency features and so on, but very few have to the science event unbinned data. We have thus excluded 

focussed on the coherence time, which is in fact a very con- the 1996 March 6th data which were obtained in a binned 

straining parameter for QPO models. In this paper, we study data mode with no spectral information. As the sensitivity 

the previously reported kilo-Hz QPOs in the X-ray burster, for QPO detection increases with the number of photons de- 

4U 1608-52, with the aim of characterizing the quality fac- tected, we have selected only those observations in which the 

tor, and hence coherence time, of its QPO over a wide range five RX TE proportional counter array (PG A) units were op- 

of frequency. We have selected this source for two main rea- erating iSradt. Rothschild, fc Swan kl 19931) . This represents 

sons: first its QPOs are known to be strong, and second in more than 250 kiloseconds of data. The dataset is divided 
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Figure 1. Time evolution of tlie 560-890 Hz lower-frequency QPO detected from 4U1608-52 over the whole data set. The image is a 
dynamical power density spectrum, in which each column is a PDS integrated over 32 s. The QPO corresponds to the boldest parts of 
the image. Note that for the longest segment of 1996 March 25th (segment 11) the QPO is barely visible in the image. 



into 21 segments recorded within different, sometimes con- 
secutive, orbital revolutions of the satellite. 

We have computed Leahy normalized Fourier power 
density spectra (PDS) between 1 and 2048 Hz over 1 s in- 
tervals (1 Hz resolution), using X-ray photons of energy be- 
tween 2 and 40 keV. It is well known that the QPO frequency 
can vary by as much as 50 Hz in a few thousand seconds. 
The goal of our analysis is to estimate Q, while keeping 
the contribution of the frequency shifts as small as possible. 
This implies that the QPO profile properties have to be es- 
timated on the shortest possible timescales. For 4U16Q8-52, 
which shows strong QPOs, this can usually be done with 32 
s. So we have averaged 32 one-second PDS and searched for 
a QPO peak using a s liding window algorithm, as described 
in lBoirin et alJ i2000l) . In segment 11, where the QPO is the 
weakest, we have averaged 128 one-second PDS. We have 
preferred to use 32 summed 1 Hz spectra rather than a sin- 
gle power spectrum obtained over, say, 32 s as in this way 
we can make a direct comparison with the results discussed 
below (§2.2), which rely on having 1 s time resolution. Using 



32 s spectra would formally give better frequency resolution, 
but we have found that the reduction in the measured widths 
and the improvement in the uncertainties in the parameters 
are very small. 

To obtain reliable fits, only features of significance 
larger than 6a are considered, although decreasing this 
threshold to 5a does not affect the results. In our data set, 
QPOs are detected on 1996 March 3rd and between 1998 
March 24th and 29th (see Figure [Hand Table [T||. 

Over the short integration times considered here, the ef- 
fects of statistical fluctuations on the observed QPO profile 
can be very important (e.g., all the power can be concen- 
trated into a single frequency bin), hence fitting of the QPO 
has to be done with great care to minimize potential bi- 
ases. The QPO is fitted over a 100 Hz frequency range (50 
Hz on each side of the QPO peak) with a Lorentzian hav- 
ing three parameters (amplitude, frequency and Av) plus 
a constant representing the counting-statistics noise level 
(generally very close to the theoretical value of 2). We note 
that the noise level can also be estimated by averaging the 



On the high coherence of kilo-Hz Quasi- Periodic Oscillations 3 



Table 1. Mean Q factor of the low-frequency kilo-Hz QPO observed from 4U1608-52 in the 21 segments of data. The count rate is given 
in cts s~^ in the 2-40 keV range for the 5 PCA units. Nint is the number of Is PDS averaged. Ns is the number of averaged PDS. Nd is 
the number of PDS in which the QPO was detected above the 6 a threshold, i^min max are the minimum and maximum QPO frequencies 
measured in the 32 s PDS. The RMS is the root mean square amplitude given as a percentage of the source count rate. Au is the fitted 
FWHM of the QPO obtained, by shifting and adding the 32 s (128 s for segment 11) PDS to a reference frequency given as fqpo. Q is 
obtained by dividing the latter by Au. 



Segment 


Date 


cts s 


Nint 


N. 


Nd 




,max 


RMS (%) 


i^qpo 


Av 


Q 


1 


03/03/96- 


-19:18 


3165.4 


32 


106 


81 


846.0- 


-891.5 


6.6±0.1 


866.6±0.1 


5.8±0.2 


150.4±5.0 


2 


03/03/90- 


-20:54 


3041.7 


32 


115 


115 


821.7- 


-840.0 


7.7±0.1 


829.5±0.0 


4.1±0.1 


200.4±4.6 


3 


CY\ /no jQ.a 


-zz.ou 


3079.0 


32 


57 


57 


827.5- 


-842.3 


7.7±0.1 


834.5±0.1 


4.2±0.1 


199.7±6.3 


4 


03/24/98- 


-15:37 


1971.9 


32 


86 


86 


749.0-773.6 


9.2±0.1 


761.5±0.0 


4.2±0.1 


181.4±5.1 


5 


03/24/98- 


-17:07 


2060.3 


32 


111 


110 


795.5- 


-820.5 


8.6±0.1 


808.6±0.0 


3.8±0.1 


212.7±5.6 


6 


03/24/98- 


-18:43 


1959.1 


32 


111 


111 


732.2- 


-779.2 


8.9±0.1 


761.5±0.0 


4.0±0.1 


192.7±5.0 


7 


03/24/98- 


-20:19 


1830.4 


32 


103 


103 


662.5- 


-710.7 


9.2±0.1 


688.5±0.1 


4.9±0.1 


139.2±4.0 


8 


03/24/98- 


-22:02 


1767.9 


32 


98 


86 


642.7- 


-669.6 


8.7±0.1 


653.4±0.1 


6.0±0.2 


108.9±3.9 


9 


03/24/98- 


-23:54 


1425.7 


32 


68 


46 


640.2- 


-667.1 


9.0±0.2 


657.4±0.1 


5.7±0.3 


116.3±5.9 


10 


03/25/98- 


-13:55 


1564.3 


32 


83 


20 


591.7- 


-609.0 


7.8±0.2 


600.5±0.2 


9.5±0.6 


63.5±4.0 


11 


03/25/98- 


-15:31 


1543.4 


128 


27 


8 


560.9- 


-580.5 


6.5±0.2 


570.1±0.6 


16.3±1.7 


35.0±3.6 


12 


03/25/98- 


-18:43 


1529.1 


32 


61 


23 


597.1- 


-633.2 


8.3±0.2 


611.1±0.2 


9.0±0.6 


68.0±4.8 


13 


03/26/98- 


-13:55 


1716.8 


32 


85 


84 


754.7- 


-788.2 


8.9±0.2 


773.5±0.1 


4.3±0.2 


178.6±6.3 


14 


03/26/98- 


-15:31 


1586.5 


32 


109 


109 


709.5- 


-793.4 


9.7±0.1 


749.5±0.0 


4.2±0.1 


176.9±5.0 


15 


03/26/98- 


-17:07 


1655.6 


32 


55 


54 


776.3- 


-799.1 


9.3±0.2 


791.4±0.1 


4.3±0.2 


182.2±7.4 


16 


03/27/98- 


-12:28 


1082.0 


32 


56 


55 


697.3- 


-745.5 


12.0±0.2 


725.5±0.1 


4.7±0.2 


153.2±6.5 


17 


03/27/98- 


-13:54 


1066.3 


32 


82 


82 


720.6- 


-752.6 


11.6±0.2 


738.5±0.1 


4.2±0.1 


177.5±6.2 


18 


03/27/98- 


-15:30 


1132.9 


32 


111 


97 


751.0- 


-797.0 


10.5±0.2 


773.4±0.1 


4.9±0.2 


159.4±5.5 


19 


03/28/98- 


-14:03 


942.0 


32 


73 


55 


767.6- 


-813.2 


ll.OitO.3 


786.6±0.1 


4.3±0.2 


183.0±8.4 


20 


03/29/98- 


-09:06 


766.3 


32 


51 


31 


664.7- 


-688.9 


12.2±0.4 


674.6±0.1 


5.8±0.4 


117.2±7.4 


21 


03/29/98- 


-10:42 


759.0 


32 


65 


41 


661.8- 


-700.0 


12.1±0.3 


679.6±0.1 


4.8±0.3 


141.6±8.0 



power in a region where no signal is expected, e.g. above 
1400 Hz. Although the two methods give very similar re- 
sults, we find that there is a trend for the latter method to 
give slightly lower A//. This is because the fitted mean level 
in the region of the QPO is generally very marginally lower 
than that above 1400 Hz^ prob ably due to subtle dead time 
effects fe.g. lZhang et al.llll995l) ). As we are interested in this 
paper in setting upper bounds to the QPO width, we have 
preferred to fit the constant level together with the QPO. 

In a more dramatic way, the results of the fits are sen- 
sitive to the way the error bars are set on the measurements 
of power. In a Leahy normalised PDS, the theoretical error 
bars on the power are P/\/N , where N is the number of PDS 
averaged. However simply using the power P in each bin to 
estimate the error on that bin leads to a bias caused by sta- 
tistical fiuctuations, because bins which happen to have low 
values of P are given high weight. We have therefore used an 
iterative procedure in which a first fit is done with uniform 
error bars and subsequent iterations are made in which the 
P used in estimating the error bars is based on the fitted 
model. The uncertainties in the fitted parameters are com- 
puted such that = Xmin + li with the other parameters 
refitted where Xmin is the corresponding to the best-fit 
value. 

2.1 Mean QPO width 

First we wish to estimate the mean width (A;/) of the QPO 
for each of the 21 segments of data, based on the measure- 
ments made every 32 s. For this purpose, we apply the shift- 
and-add technique (e.g., Mendez et al. 1998) to the 32 s 
PDS. Each 32 s PDS is fitted to obtain the QPO frequency. 



then shifted to a reference frequency. The resulting PDS is 
then averaged to produce a single PDS, from which the mean 
QPO width (A;/) can be fitted. The reference frequency cho- 
sen is the mean QPO frequency for the segment of data 

In those 32 s PDS where the QPO is not detected above 
6(T, the QPO frequency used for the shifting of the PDS 
has been estimated by linearly interpolating between neigh- 
bouring detected frequencies on either side. Similar results 
are obtained if we use a spline interpolation. We have also 
checked that increasing the PDS averaging time (say to 64 
seconds) to increase the number of QPO detections yielded 
consistent results within error bars for A;/. 

This method should overestimate somewhat the true 
because the frequencies determined every 32 s and used 
in the shift-and-add have typical error bars between 0.5 and 
1 Hz, introducing a blurring of the QPO profile. Obviously 
any variation faster than 32 s would also broaden the QPO. 
The mean quality factor (Q), obtained by dividing i/qpo by 
Au should thus be considered as a lower limit. Nevertheless, 
we find values of Q larger than ~ 150 in 14 out of the 21 
segments analysed, and values as high as 213 are found, 
indicating that the large Q factor reported by Berger et al. 
(1996) is by no means exceptional. The results are listed in 
Table 1. 



2.2 Attempting to correct for the observed drift 

Even with PDS averaged over only 32 s, the measured width 
may still be affected by the frequency drift. As can be seen 
in Figure |5| drifts as large as 5 Hz between consecutive 32 
s PDS are not unusual. Can we estimate what would be the 
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Figure 2. Time evolution of the QPO frequency in segment 2 as 
followed every 32 seconds. Frequency drifts of up to 5 Hz in 32 s 
are observed. 



Figure 3. Time evolution of the QPO frequency in segment 5 as 
followed every 1 second using a sliding window of 32 s. 



width of the QPO in a one-second PDS (i.e. with negligible 
contribution from the frequency drift)? It is impossible to 
fit directly the one-second PDS but one can still make some 
first order estimates. 

If one assumes that the QPO frequency drifts linearly, 
one can remove the drift contribution using the quadratic 
relation 

Au^ ~ AuiB^ + \udrm\ (1) 

as an approximation, where Au is the mean QPO measured 
on 32 s with the method described above, Auis is the mean 
width expected in a 1 s PDS and li/driftl is the mean absolute 
frequency differences between consecutive 32 s PDS. 

Alternatively, one can try to account for the drift by 
again using the shift-and-add technique, applied this time 
to the 1 s PDS. For this, we must estimate the QPO fre- 
quency every second. To do so, here we have used a sliding- 
window technique; the frequency is estimated by fitting the 
QPO from a PDS which is the average of 16 one-second 
PDS before and 16 one-second PDS after the time bin (to 
avoid any possible biases, the PDS of the central time bin 
is excluded in the average). Then each 1 s PDS is shifted 
in frequency and averaged to produce one single PDS from 
which the QPO can be fitted to obtain Ai^sa,is. We have 
checked through simulations of a pulsar-like signal, drifting 
in frequency by about 4 Hz per 32 s (to reproduce the width 
measured in segment 2, Table0, that the method described 
above is accurate and able to recover the pulsar nature of 
the signal (in the shifted PDS, the power is concentrated 
into a single frequency bin). 

For maximum accuracy, we have applied this method 
to the segments of data in which the QPO is almost always 
(more than 90% of the time) detected in 32 s time intervals 
(segments 2-6, 13-18). An example of the time evolution of 
the QPO frequency measured every second using the sliding 
window technique is shown in Figure The results of the 
fits to the shifted and added 1 s PDS are listed in Table 121 
together with |i^drift l and Auis from equation As can be 



Table 2. The results of shifting and adding the PDS on a 1 s 
timescale for those segments of data where the QPO is detected 
in more than 90% of the 32 s time intervals. Ai/=Ai/sa.32sis the 
same as Table 1. |i^driftl is the mean absolute frequency drift be- 
tween consecutive 32 s averaged PDS, Auis is the width expected 
assuming a linear drift within the 32 s time interval, and At'sa.is 
is the fitted width of the QPO in the shifted and added PDS. 



Seg 


Aj^=Al/sa,32s 


kdrift 




A!^sa,ls 


2 


4.1± 0.1 


1.3 


3.9 


4.2± 0.1 


3 


4.2± 0.1 


1.2 


4.0 


4.3± 0.1 


4 


4.2± 0.1 


1.4 


4.0 


4.2± 0.1 


5 


3.8± 0.1 


1.4 


3.5 


3.9± 0.1 


6 


4.0± 0.1 


1.4 


3.7 


4.0± 0.1 


7 


4.9± 0.1 


1.7 


4.7 


5.2± 0.1 


13 


4.3± 0.2 


1.5 


4.1 


4.7± 0.2 


14 


4.2± 0.1 


1.6 


3.9 


4.5± 0.1 


15 


4.3± 0.2 


1.6 


4.0 


4.4± 0.2 


16 


4.7± 0.2 


2.2 


4.2 


5.2± 0.2 


17 


4.2± 0.1 


1.9 


3.7 


4.4± 0.2 



seen, the recovered Au are systematically larger than would 
be expected from equation Q and are actually larger than 
or equal to those obtained by shifting and adding the 32 s 
PDS. 

There are two possible effects that may combine to ex- 
plain this quite surprising result. First, the individual fre- 
quency estimates have typical error bars (between 0.5 and 
1 Hz) of the order of the gain expected in the shift-and-add 
process. These errors will naturally introduce a blurring of 
the QPO profile in the shifted PDS. We have checked this 
effect using simulations of a QPO signal of known frequency, 
and similar amplitude and width as the real data. Alterna- 
tively, this might be an indication that the phase, ampli- 
tude, or frequency of the underlying QPO signal vary on 
time scales significantly shorter than 32 s. 
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Figure 4. Top) QPO coherence time as a function of QPO RMS 
for segment 1, when the QPO disappeared at about 890 Hz. (Bot- 
tom) QPO Coherence time versus QPO frequency for the same 
segment of data. To reduce scattering in the data points, the av- 
eraging time for the PDS is 64 seconds instead of 32 s. 



2.3 Coherence time of the oscillator 

In this paper we define the coherence time of the oscilla- 
tor on the assumption that the signal consists of sine wave 
shots with exponentially decaying amplitudes with a time 
constant r. Such a signal will produce a Lorentzian with a 
FWHM Au^1/{ttt) in the PDS. 

Berger et al. (1996) reported a correlation between 
Aiy and the frequency, as well as an anticorrelation between 
the RMS and frequency, using data from the first segment 
analysed here and a PDS integration time of 100 seconds. As 
can be seen in Figure this segment of data is remarkable 
because the QPO reaches its highest frequency and seems to 
weaken as it does so. In Figure 2] we show the dependence 
of the inferred coherence time versus QPO RMS and QPO 
frequency. There is a strong anti-correlation between the 
frequency and coherence time (Spearman correlation coeffi- 
cient of -0.75 corresponding to a null-hypothesis probability 
of 1.9 10~^), and a weaker correlation between the coherence 
time and RMS (Spearman correlation coefficient of 0.28, cor- 
responding to a null-hypothesis probability of 0.06). 

We have investigated whether these trends were present 
in the whole data set, in which the frequency spans a much 
wider range. The relationship between the coherence time 
and fqpo for all the data analysed here is shown in Figure 
|5] To compute the coherence time, we have used Av from 
Table El There is a clear positive correlation between the 
two quantities until j/qpo reaches a maximum level around 
~ 800 Hz (albeit with some scatter). After the maximum, an 
anticorrelation between the QPO frequency and coherence 
time, of the type seen within segment 1, seems to be present. 
A similar behavior was reported from 4U1636-53 for its lower 
kilo-Hz QPO (di Salvo et al. 2003; see also van der Klis et 
al. (1997) for the upper kilo-Hz QPO of Sco X-1). 

The overall relationship between the RMS amplitude of 
the QPO and its coherence time is shown in Figure |S| The 
behaviour is much more complex; both positive and negative 



Figure 5. Dependence of the QPO coherence time on frequency. 
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Figure 6. Dependence of the QPO coherence time on RMS. The 
numbering refers to Table 111 



anticorrelations are seen. Again, the pattern seen in segment 
1 is not typical. 



2.4 Coherence of the higher-frequency QPO 

We have searched for the second QPO in the data using the 
shifted 32 s PDS used above. We detect an otherwise in- 
visible higher-frequency QPO in 5 of the 21 segments used 
here. A weaker signal is detected in segments 1 to 3, 13 to 
15, 16 to 18 and 20-21, which we group to obtain a mean- 
ingful fit. The results of fitting the higher-frequency QPOs 
are listed in Table |3 Given the relatively large width of the 
high-frequency QPO, the window for the fit is 200 Hz on 
each side of the QPO peak, excluding the region around the 
lower-fre quency QPO. These Q P Os hav e been previously re- 
ported in lMendez et al.l ljl998allblll999l) . What is noticeable 
is that the higher-frequency QPO is much broader than the 
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Table 3. The higher-frequency QPOs detected after applying the shift-and-add procedure to the 32 s PDS. i^ahift is the frequency to 
which the 32 s PDS were shifted. Ai^io„ is the width of the low-frequency QPO. I'high, At^igh and RMShigh are respectively the frequency, 
the width and the RMS of the high-frequency QPO. To obtain reliable fits, segments 1 to 3, 13 to 15, 16 to 18 and 20-21 were grouped 
together. The errors are again computed as = X^in ^^ 



Seg J^shift At-iow i^high Ai^high RMShigh (%) 



1-3 


843.5l°:0^3 


^"^-0.1 


1103.9+^7'^ 


179 e+^'i-^-'^ 


4.1±0.2 


7 


688.5lH^ 


4 g+O-i 


994.0t^;J ■ 


26.1+9^^7^ 


3.4±0.7 


8 


653.41^^8 


D.U_g 2 






2.9±0.8 


10 




9 5+°-^ 


9OI.8I7 




6.4±0.5 


11 


57Q.lt°:ll 


16.3+1-^ 


879.51^;^ 


89.8l?°:« 


8.4±0.4 


12 


611 l+°-2^ 


9 0+°-^ 
^■'-'-0.6 


930.7+^7 
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lower one, with Q less than ^ 10. T his is a well known ob- 
servational fact (Ivan der KlislEoOoll . Interestingly, there is 
no apparent correlation between the coherence time of the 
higher frequency QPOs and its frequency (of the type seen 
in Figure 13 for the lower-frequency QPO), but we note that 
the frequencies do not span a wide range and that the error 
bars are relatively large. 



3 DISCUSSION 

We have carried out a systematic study of the kilo-Hz QPO 
of 4U1608-52, with the aim of investigating the coherence 
time of the underlying oscillator. By following the lower kilo- 
Hz QPO over the shortest timescales permitted by the data 
statistics (down to 32 s) to minimise the contribution of the 
long-term frequency drifts, with a shift-and-add technique, 
we have found a lower limit on the quality factor larger than 
~ 150 in 14 of the 21 segments of analysed data. We have 
shown that in this object Q is as high as 200 for thousands 
of seconds. In many sources, the strength of the signal does 
not allow the study of the QPOs on such short timescales, 
which may explain why the Q reported in the literature are 
generally lower than this. 

Frequency drifts as large as 5 Hz are observed between 
consecutive 32 s PDS. In an attempt to remove the contri- 
bution of this short-term drift to the width of the QPO, we 
have applied the shift-and-add procedure to the 1 s PDS, 
using a sliding-window technique to estimate the QPO fre- 
quency every second. If the frequency drift were linear over 
32 s, this method should remove the drift contribution to the 
width measured on 32 s and provide an estimate of the QPO 
width on a 1 s timescale. Interestingly enough, the recovered 
QPO width is not significantly smaller than the width mea- 
sured on 32 s. This surprising result could be explained by 
the fact that the error bars on the frequency are of the or- 
der of the expected gain, causing an artificial blurring of the 
QPO profile in the shift-and-add process. Alternatively, this 
might be indication that the frequency, amplitude, or phase 
of the underlying signal varies on time scales much shorter 
than 32 s. In this case, this would most likely imply that 
the Q factor reported in Table^with the present technique 
underestimates the intrinsic Q of the QPO. 

Determining the coherence time of the underlying os- 



cillator from the width of the QPO over the wide range of 
frequency spanned (between 560 Hz and 890 Hz), we have 
shown that there is a clear pattern in which the coherence 
time increases with frequency up to a maximum at ^ 800 
Hz. At both ends of the frequency range, the decreasing 
strength of the QPO seems to be associated with a decrease 
of both its amplitude and its coherence. 

We now discuss the implications of the above results 
for QPO models. There are two different classes of models 
of the kHz QPOs: those i nvolving clumps within o r abov e 
the accretion disk (e.g.. iMille r. Lam b, fc Psaltij il99d) ; 
ISteUa fc Vietril lll999lllvan der Khj l|2004^ for a recent re- 
vi ew), and those i nvolving oscillations of the disk (reviewed 
in lWagoneJ ^199^ . and .Katg 1,2001)). 

3.1 Clumps 

If inhomogeneities can form in the accretion flow in the form 
of clumps which are more luminous than their surroundings 
and which orbit around the central star, the viewing geom- 
etry might be such that the clumps will produce luminosity 
variations (.Bathiiil973i) . Inhomogeneities lying at a radius R 
and having a radial extent AR will be sheared by differential 
rotation and will therefore have a finite lifetime. The maxi- 
mum lifetime for such a clump to be stretched to an axisym- 
metric ring is given by Ts = (2/S )(R/ AR)Pk, where Pk is the 
Kepl erian period at radius R iBath. Evans, fc Papaloizoul 
I1974I) . In the clump model, this lifetime is comparable to 
the coherence time of the QPO r ~ r^. Such clumps will 
produce a lum inosity variation of at most AL/L = AR/ R 
l|Pringlell98'll) . Hence, Ts cx R/AR oc (AL/L)"\ so high co- 
herence of the signal and its high amplitude are incompatible 
in the clump model. In our case AL/L is of the order of 10% 
(see Figure [SJ. This would imply ~ lOPfc < 0.01 second 
(or Q ~ 30), much shorter than the coherence time we found 
in our analysis (see Figure (SJ. It is also worth noting that 
in the above formula the clump lifetime (and signal ampli- 
tude) are predicted to decrease with increasing frequency, 
(e.g. Livio fc Bath 1982), unlike the observed behaviour in 
Figures |21 and |S| 

This is clearly not what we observe in this source 
(see Figure |5J|. Alternatively, one can make the assump- 
tion that the radial extent the blob is of the order of 
the disk height {h). Shakura and Sunyaev (1973) estimate 
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/i ~ 3 X lO''cm(L/LEdd)(-M"/M0), and the radius of the rele- 
vant orbit is on the order of the stellar radius, or a few times 
larger. For an accreting neutron star with a large fraction of 
Eddington luminosity, h/R ~ 1/10 is a reasonable estimate, 
which again leads to the same conclusions. 

We have also found that the higher-frequency kilo-Hz 
QPO is much fainter and is characterised by a much lower 
quality factor. This result is hard to reconcile with models in 
which the higher- frequency QPO is a Keplerian frequency at 
some radius (e.g., the sonic point. iMiller. Lamb, fc Psaltij 
il998l) l. and the lower QPO is a beat frequency generated 
by interaction of the neutron star radiation with clumps 
originating at that particular radius. This model predicts 
that the width of the two QPOs should be similar, because 
the stellar spin frequency is expected to be highly coher- 
ent. In 4U1608-52, we have shown the width of the higher- 
frequency QPO to be much larger than the width of the 
lower- frequency QPO. 

Another model of QPOs has been proposed in which 
clumps leave the accretion disk and for some reason follow 
test particle orbits. In this so-called relativistic precession 
model (e.g. Stella & Vietri 1999), the clump is responsible 
for all three frequencies: two kilo-Hz QPOs and a lower- 
frequency QPO, also called the horizontal branch oscillation. 
Since the same agent is responsible for the two kHz QPOs, 
the coherence times should be about the same. We have 
shown here that the coherence time of the higher-frequency 
QPO is at least a factor of ten lower than that of the lower- 
frequency QPO (See Table 3). 

More stable non-linear structures (vortices) may persist 
in the accretion disk, but these will be subject to radial 
inflow unless they are anchored by an external magnetic fiel d 
jAbramowicz et alll993lLovelace et al.ll99gl : lTaggeil200ll) . 
The inflow velocity at radius r in an a-disk of thickness h, 
is a factor of a{h/r)^ lower than the orbital velocity, where 
Q > 10~^ is the dimensionless viscosity. So, for typical disks, 
within ~ 10 orbits the fluid is carried inwards by a distance 
equal to a few per cent of its initial radius, and suffers a 
percentage change in orbital frequency which is 1.5 times 
larger. 

Thus all models of this class seem to be incompatible 
with a high Q factor, and the relationships found between 
the coherence time with frequency and RMS amplitudes. 



3.2 Disk oscillations 

There is growing evidence that the oscillations originate in 
the disk itself, in particular because the same phenomena 
are observed in a wide range of a ccreting X-ray sources, 
from white dwarfs to black holes jMauchdl200a) . Several 
models, relating QPOs to disk oscillations have been pro- 
posed, some of them making clear predictions of what the 
quality factor of the QPO should be. This is the case, 
for instance, for the transition l ayer m odel proposed by 
iTitarchuk. Lapidus. fc Muslimovl il998h . In this model, a 
shock occurs where the Keplerian disk adjusts to the sub- 
Keplerian flow and the transition can undergo various types 
of oscillations under the influence of the gas, radiation, mag- 
netic pressure, and gravitational force. In this case Q can in 
principle reach 100. It is e xpected to increase with frequency 
but also with luminosity iTitarchuk. Lapidus. fc Muslimovl 



il998ll l. In our data set it is clear that there is not a good 
correlation between Q and the source count rate. 

For acoustic wave s in isotropic turbulence 
iGoldreich fc Kumad lll98d) predict Q ~ h/{Xa). Here 
a ~ 0.1 is viscosity parameter for a Shakura-Sunyaev disk 
and the wavelength. A, is much greater than h, yielding 
a rather low value of Q. Anisotropic turbulence increases 
this to Q = X/{ha) ~ 10, as observed in black holes. For 
g-modes and isotropic turbulence, Q is the higher value of 
1/q and y^r/h, but the highest values of Q are obtained 
for g-modes in anisotropic turbulence, {r/h)/a ~ 100. For 
c-modes, which correspond to a warp revolving at the Lense- 
Thirring frequency (a general-relativistic effect), this value 
is enhanced by a factor approximately equal to the square 
of the ratio of orbital velocity to the speed of sound (the dis- 
cussion in this section follows lOrtega-Rodn'guez fc Waeoneil 
(2000)). In the very hot disk around the neutron stars this 
enhancement may still be insufficient to explain the high 
Q reported here. However, some modes are excited, rather 
than d amped, by viscosity and may enter the non-linear 
regime dOrtega-Rodrfguez fc Wagoneillioool) . 



4 CONCLUSIONS 

We have shown that the kilo-Hz QPO in 4U1608-52 is a 
highly coherent signal with an inferred coherence time of 
up to about one tenth of a second. This result, as well as 
the dependence of the QPO coherence time on frequency 
and amplitude which we report are hard to reconcile with 
QPO models involving orbiting clumps. This leaves the ac- 
cretion disk oscillations as the the most likely alternative for 
accounting for these signals. 
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